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ABSTRACT 


Electricity constitutes the input into many products that produced by industry and used by people. Hence, 
it can be considered as a product or service that has vital importance in human life and economy. Since 
it has such special properties of instantaneous production and consumption obligation and unfeasible 
storage, electricity market is not like other markets. In a competitive electricity market, generation com- 
pany faces price risks and delivery risks. So that risk management is an important part of a generation 
company and can deeply effect companies’ profitability. This paper focuses on electricity generation asset 
allocation between bilateral contracts, such as forward contracts, and daily spot market, considering con- 
straints of generating units and spot price risks. The problem is to find the optimal portfolio based on 
known electricity generation total costs, bilateral contract prices, it employed Turkish historical balanced 
market hourly system marginal and day-ahead hourly market prices between of 2006 and 2011. There 
are limited studies about portfolio optimization in electricity markets in literature and this paper should 
be considered frontier study taking spot market’s hourly prices separately as risky assets. Markowitz 
mean-variance optimization which is claimed to be the beginning of modern portfolio theory in financial 
sector is used to demonstrate this approach. Mean-variance optimization has been successfully applied 
to all cases that modeled for electricity market. Some suggestions for future work are also listed in this 
paper. 

© 2011 Elsevier Ltd. All rights reserved. 


Mean-variance optimization 
Risk management 
Turkish electricity market 


Turkey 
Contents 
V3. INMOdMUCHON: sccss:is ccc aids noan cate ssiesbeadinentdas deste cents weed et eeawead E NaN Ese casGe DbeDaNORe daa deeU TED CERES wa ew ee eae 
2: Current:status of Turkish electricity Market iicc..cccicsscctens seekers idek kaai weds atre EN E a KEN AGa ADEE E dead aes swatare geet 
3. Portfolio optimization theory srece caters sate ovsdcsasiedetenaccadains otaaes E A E E ER R sw'ee doses ceateepeens dem RE 
331... Modern portfolio theory. ieceri jeehacae sencisis cassie Load EErEE E RE IET TES E E EA NE A tale E ERS 
3.2. Portfolio selection and determination of utility function. .............ssessssesresesesrrsserrsresesresrsresrsrerreresrssrsreresrererresesre: 
4. Data andimethodolosy::.. cassadcctsiis to aaea nE a E EE va catia EEA E E a aa ON EEE EE E E E 


4.1. Determination of risky and risk-free assets .......... o...n eeneneenesesosresorrereresreseseresseresosreresrererrerereesesceresesersresesrerre: 
4.2. Construction of three base MOdelS............c:ceceeennenceceuecceccccccscesccennenceennenesetseecenbbeneetenenenteneeensnbeseuteneceecesens 
4.2.1. Standard mean-variance optimization model .............. cc cece cece eee eee eee e ene e eee e ee ee eens eee beeen eee eeeeeeeeeenees 
4.2.2. Mean-variance optimization model with upper investment constraints 
4.2.3. Mean-variance optimization with one risk-free asset ..............ccceeeeeeeee eee e cece eens 


5.1. Efficient frontier and optimum portfolio for hydraulic power plants .................cceeeeeeeeeeeeee 
5.2. Efficient frontier and optimum portfolio for lignite fired thermal power plants 


* Corresponding author at: Ankara University, Faculty of Political Sciences, Department of Management, 06600 Cebeci-Ankara, Turkey. 
E-mail addresses: fgokgoz@politics.ankara.edu.tr (F. Gökgöz), meteemin.atmaca@euas.gov.tr (M.E. Atmaca). 
1 Ankara University, Nuclear Engineer & Project Manager at EUAS (Electricity Generation Company Inc. of Turkey), Ankara, Turkey. 


1364-0321/$ - see front matter © 2011 Elsevier Ltd. All rights reserved. 
doi:10.1016/j.rser.2011.06.018 


358 F. Gökgöz, M.E. Atmaca / Renewable and Sustainable Energy Reviews 16 (2012) 357-368 


5.3. Efficient frontier and optimum portfolio for natural gas combined cycle power plants..............c cece eeeee cece cece ee eeeeeeeeneneeeeees 366 


5.4. Analysis of results ...........c.cccccccesesesceseseacacacaceaseenes 
6; CONCLUSIONS. s ctrrnvctdiaseciderwecdas EDE E EEEE eoedaes E 
References): ssecsdssnddeaniiediaararveahsddeatweadas seasoned debacad deta sad 


1. Introduction 


In many countries today, there is a remarkable tendency for 
deregulation and restructuring of electric power industries. Turkey 
is one of these countries and it is also in a restructuring period. It is 
aiming to deregulate all electricity industry until 2015. However, 
this new market environment brings along some risks. In a com- 
petitive environment, generation companies’ main objective is to 
maximize their profit and minimize the associated risks, therefore 
clear determination of the risks and taking necessary steps are very 
important. It is essential to create a trading strategy with risk man- 
agement before bidding into spot markets (high profit with volatile 
prices) and bilateral contract market (moderate profit with stable 
prices). At this point, main step is the determination of portfolio 
weights and optimization of the portfolio with respect to investors’ 
risk aversion. 

Various risk management techniques have been applied to elec- 
tricity markets. Risk management can be divided into two main 
sections; risk control and risk assessment. In the study of Liu and 
Wu [1], hedging and portfolio optimization have been defined as 
risk control techniques. On the other hand risk measurement and 
asset valuation are described as risk assessment techniques. Hedg- 
ing is necessary to offset risk of position (generally in spot markets) 
and to reach this goal investors must buy some other derivatives 
like forward contracts, future contracts, options and swaps, etc. 
[2,3]. In electricity markets hedging the risk of spot price with for- 
ward contract and futures contract have been investigated [4,5]. 
Forward contracts, purchased by a consumer are commitment by 
the suppliers to provide a specified amount of electrical energy at 
a future time, are particularly useful for coordinating supply and 
demand side operations decisions under spot pricing [4]. A future 
contract obligates each party to buy or sell a specified amount 
of energy at a specified price. But, buyers and sellers of future 
contracts deal with an exchange, not with each other [6]. Hedg- 
ing techniques using futures contracts in electricity market appear 
to have lower standard deviation or risk [5]. Options, Swaps and 
special electricity contracts (callable forward contract, puttable 
forward contract) are also determined as other tools for hedging 
strategies. 

Portfolio optimization, listed among risk control techniques, 
refers to optimally allocation of energy trading tools with the aim 
of maximizing benefits while minimizing the corresponding risk. 
Except Monte Carlo Method in reference study of Vehviläinen and 
Keppo, there are two types of methods can be used to solve portfo- 
lio optimization problems: Decision analysis and modern portfolio 
theory [6,7]. Decision analysis is the determination of all possible 
events with its consequences and probabilities and the constitu- 
tion of decision tree with respect to these data and the evaluation 
of tree [8]. Modern portfolio theory (MPT), which is the subject of 
this paper, is the other technique that can be used for portfolio 
optimization [1,2,6-12]. Mean-variance optimization is the essen- 
tial part of this approach. Only few studies have been done by using 
MPT in electricity markets. It is very popular in financial literature 
but it is brand new approach for electricity markets. Risk assess- 
ment techniques for portfolio optimization are not scope of this 
study so it is not necessary to give additional information on them. 

Section 2 introduces the background of the electricity market 
with different trading regimes and pricing system which includes 
trading environment. Section 2 also introduces current status of 


Turkish electricity market. Section 3 describes the fundamental 
theory and methodology (especially mean-variance optimization) 
to portfolio optimization which can be applied to electricity mar- 
kets. Section 4 explains data and methodology and includes the 
application of mean variance to electricity market environment. 
Section 5 demonstrates results of study in graphically and table 
representation forms. Finally, Section 6 provides conclusions. 


2. Current status of Turkish electricity market 


Electricity markets in the world generally offer two types of 
market structure: spot market and physical market. In addition to 
these, there are derivative markets too. Spot markets include bal- 
ance and/or day ahead spot markets while physical markets include 
bilateral and/or physical forward contracts. There are also deriva- 
tive markets (for financial instrumentation of electricity with, 
futures, option, swap and special other derivatives) in some of the 
electricity trading regimes. 

To keep interconnected system in a stable condition, it is nec- 
essary to keep frequency value of network in a bound around a 
specific value.? Instantaneous load changes must be compensated 
by load rejection/acceptation of generation units. Excess power in 
networks causes increase in frequency while opposite case causes 
decrease in frequency. Balance markets are founded to respond this 
need in a short time (generally less than 15 min) so that they are 
defined as real time markets. As to day-ahead market all energy 
traded approximately 24h (in Turkish electricity market between 
12 and 36 h before real consumption time) before real consumption 
time. The common scenario among these markets is that they all 
involve a centralized auction mechanism, by Independent System 
Operators (ISOs), Regional Transmission Operators (RTO) or any 
such organizations, to determine which generation units should be 
deployed and how much energy should be produced to meet the 
demand by selected unit [1]. Depends on the countries spot mar- 
ket mechanism preferences, trading intervals to submit bids in spot 
market can be different. For instance, in NEM of Australia, the time 
interval for electricity trading is defined as 30 min but all counter- 
parts can give their bids for every 5 min. In Greece, STA has single 
and separate market for each hour. In England and New Zealand, 
NETA and NZEM have 48 half hour trading periods for each day. 
In IMO of Ontario bids and/or offers are submitted for every 5 min. 
As to Elspot in Nordic market, market participants bid for purchase 
and sale of power contracts 1-h duration for all 24h of the next 
day. Turkish electricity market has been operated for hourly trading 
interval, too. 

Three types of pricing mechanism are common in world’s energy 
markets: Uniform marginal pricing, zonal pricing and locational 
(or nodal) pricing. In balance market, some other pricing system 
which is called as pay-as-bid or discriminatory pricing system can 
also be seen for instance Wales and England market. In markets, 
where uniform marginal pricing were adopted, only one energy 
price is used for ex post settlement for each trading interval. 
Generation companies can make certain of their revenues by sign- 
ing bilateral contracts with their customers at fixed energy prices. 
In zonal pricing systems, there is a geographical point of view to 


2 American network system is being operated at 60 Hz frequency and 110V, Euro- 
pean and Turkish network systems are being operated at 50 Hz frequency and 220V. 
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Fig. 1. European market model/multi environment advanced market model. 


market structure and all system divided sub systems with respect 
to pre-defined geographical areas, called zones. When there is no 
congestion between zones, one uniform market clearing price is 
used throughout the system. Under congestion conditions all zones 
separated each other and produce their own market clearing prices, 
called zonal prices. These pricing systems have congestion risks for 
bilateral markets because of their uncertain and unpredictable pric- 
ing mechanisms. Zonal pricing is still used in Nordic (11 zones) and 
some other markets [13]. Locational marginal pricing is the deter- 
mination of prices in each location or node of power system. When 
there is no congestion, one market clearing price is used in the sys- 
tem (ignoring transmission loses) same as zonal pricing. In zonal 
and locational pricing regimes, except contract with customers at 
the same location or zone, all contracts are risky [1]. Theoretical 
point of view, methodology for each pricing system is the same. 

In Turkey, deregulation and reconstruction processes are pro- 
ceeding. Turkish authorities are planning to reach completely 
deregulated market environment before 2015. With respect to 
Turkey Progress Report of European Union 2009, the thresh- 
old for eligible consumers was further reduced to 480,000 kWh, 
equivalent to opening 50% of market [14]. After this report Turk- 
ish Energy Market Regulator Authority published new threshold 
value as 100,000 kWh [15]. Turkey is adopting European Market 
Model to itself, in that respect Nordic Electricity Market struc- 
ture is a good example for Turkey. Turkish electricity market 
consists of two main market environments now: balanced mar- 
ket for balancing load imbalances as a spot market and bilateral 
contract application. Day-ahead planning is in service to pre- 
pare market participants to day-ahead market. In near future (at 
the end of 2011) day-ahead market will be available for energy 
traders. In Fig. 1 common details of target model of Turkey can be 
seen. 


3. Portfolio optimization theory 
3.1. Modern portfolio theory 


In classical portfolio theory, it is thought that you can 
avoid your investment risk by constructing your portfolio from 
different industries’ shares, treasury bills with different maturities, 
and different foreign currencies, in other words via diversification. 


In modern portfolio theory, taking into consideration of relations 
(correlations, co-movement) of portfolio assets with each other is 
very important and diversification is necessary but not the only 
solution of portfolio risk management. Taking these co-movements 
of securities into account satisfies an ability to construct a port- 
folio that has the same expected return and less risky portfolio 
than a portfolio constructed by ignoring the interactions between 
securities [1]. 

It is impossible for investors to know exactly what the asset’s 
return and risk will be tomorrow, so they are using some forecast- 
ing tools to determine these values. The main problem by which 
investors must be faced is to figuring out weight percentages of 
assets being in the portfolio. This problem is called as portfolio 
selection problem. To find a solution to this problem, Harry M. 
Markowitz published a paper, which is called as milestone study 
for portfolio theory and fundamentals of modern portfolio theory, 
“Portfolio Selection” in 1952 [3]. Markowitz, in his famous paper, 
argued that a process of portfolio selection can be divided into two 
stages. “The first stage starts with observation and experience and end 
with beliefs about the future performances of available securities. The 
second stage starts with the relevant beliefs about future performances 
and ends with the choice of portfolio.” His paper is concerned about 
the second stage [11]. His body of work on portfolio theory resulted 
in him being awarded Nobel prize in economics in 1990. He shared 
this prize with Merton H. Miller and William F. Sharpe and they 
were awarded jointly for their pioneering work in the theory of 
financial economics. After Markowitz’s famous paper, theory was 
amplified by Sharpe in 1964 and by Linther in 1965 [16,17]. The 
addition of a risk-free asset by Sharpe and Lintner in the mid 1960s 
led to the capital market line and the CAPM [18]. 


3.2. Portfolio selection and determination of utility function 


Markowitz’s portfolio theory is based on a mean-variance opti- 
mization process that searches for efficient portfolios. An efficient 
portfolio means one that provides minimum risk for a given level of 
return or maximum return for a given level of risk [19]. It must be 
done another basic assumption that rational risk-averse investors 
select their portfolios using only mean-variance criteria. The main 
assumptions of the mean-variance analysis are based on the fol- 
lowing issues [20]: 
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Fig. 2. Efficient frontier and efficient portfolios. 


e All investors are risk averse so that they prefer less risk to more 
for the same level of expected return. 

e Investors have the information regarding the expected returns, 
variances and covariances of all assets. 

e Investors need only to know the expected returns, variances and 
the covariances of returns to determine optimal portfolios. 

e And there exist no transaction costs or taxes limitation. 


The required inputs necessary for mean-variance optimization 
model are the expected returns for each asset, the variance of each 
asset, and the covariances between assets. Ifthe asset’s returns obey 
anormal distribution, then the entire distribution of a portfolio can 
be described by the mean and variance only [21]. Markowitz mean- 
variance optimization process produces an efficient frontier (Fig. 2) 
which consists of efficient portfolios sets on it. 

The expected return for a particular portfolio, which includes 
“N” assets, can be expressed as: 


N 
E(tp) = X Xiri (1) 
i=1 


while “N” describes number of assets that take place in portfolio, 
“X” denotes the proportion (weight percentage) of ith asset in port- 
folio and “r;” denotes expected return of ith asset. The variance for 
a particular can be expressed as: 


OF = DD Xie (2) 


“oj” is the covariance between the returns on the ith asset and the 
jth asset. The covariance oj; measures how many the returns on 
two assets move together. 

The basic mean-variance optimization model includes the min- 
imization of the portfolio’s variance under three fundamental 
constraints. These constraints are expected return of portfolio 
must be equal to target return; the sum of the proportions of 
financial assets that are being in portfolio must be equal to “1” 
and finally non-negativity condition for assets’ proportions. Under 
above determined conditions it can be set the basic equation as: 


N N 
Min. (ož) = SOS XXoi (3) 


i=1 j=1 


s.t. 


N 
X Xr =Te (4) 
i=1 

N 
Sox; =1 (5) 
i=1 


Xi>0, YXieli=1,2,...,N] 


After solution of this problem set it can be reached the efficient 
portfolios and efficient frontier as shown in Fig. 2 for N risky assets 
portfolio. 

The main issue for the basic mean-variance optimization model 
is the determination of optimal proportional allocation “X;” to the 
ith financial asset. After determination of efficient portfolios and 
efficient frontier, it is important the determination of investor’s 
utility functions that represent investor’s risk aversion level. It is 
assumed that each investor can assign a utility score to compet- 
ing investment portfolios based on expected return and risk of 
those portfolios. Combining two, we can define the objective or 
utility function for an investor in terms of expected return E(r) and 
variance of returns g? as follows [1,2,6,22,23]: 


U=E(r)— sho? (6) 


where U is the utility value and “A” is an index of investor’s risk 
aversion. The factor of 1⁄2 is just a scaling convention. In financial 
approximations, A=3 is taken for representation of average risk 
aversion, A>3 means more risk aversion and A<3 means less risk 
aversion. In the light of above information, now it can be found 
optimum portfolio for an investor that has level “A” risk aversion, 
after the solution of below equations set: 


U =E(rp) — 5403 (7) 


X;>0, VX;eli=1,2,...,N] 


where 


4. Data and methodology 
4.1. Determination of risky and risk-free assets 


To apply mean-variance optimization to electricity market, first 
of all it must be determined risky assets and to enlarge this appli- 
cation also we need a risk-free asset. In Turkish electricity market, 
we have system marginal prices and day-ahead market prices since 
August 2006. By knowing electricity generation costs with respect 
to their types (Hydraulic Power Plant, Lignite fired Thermal Power 
Plant, Natural Gas Combined Cycle Power Plant), we will use these 
data to represent risky assets. Bilateral contract prices with fixed 
values (under the guarantee of clearing authority) will be used to 
represent risk-free asset. Turkish electricity market's daily average 
spot prices frequencies between August 2006 and April 2011 can be 
seen in Fig. 3. We assumed that price frequency has a normal dis- 
tribution around mean. In calculations, yearly average (US Dollar 
selling) currency prices have been used [24]. 

Spot market electricity prices and the amount of electricity 
consumptions (demand for electricity) have volatile values but 
they move together with high correlation. Both of these two 
index can importantly change with respect to different hours of 
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Fig. 3. System marginal prices and day-ahead prices frequency distribution 
between August 2006 and April 2011. 


day, weekdays, weekend, months, seasons, years, and also change 
supply-demand balance and system congestions but each of them 
follows a characteristic way. The most well-known characteristic 
is hourly demand characteristics of consumers. Starting from this 
point that consumers have similar consumption characteristics in 
the same hours of a day, by using spot prices of 24h of day, 24 risky 
assets have been determined. Above mentioned assumptions have 
been taken to support this approach: 

Daily load curves that have chosen randomly for different years, 
seasons, months, and weekdays can be seen from Fig. 4. In spite of 
the all differences, consumers have similar consumption profiles for 
five different sample days. In such a way that load curves generally 
decreases after 10 p.m. and increases after 7 a.m. The amount of 
electricity consumption is low for the first 8 h of a day. In all curves 
it can be seen two or three head/shoulder formation shapes. Except 
extreme conditions like public holidays, religious holidays, very hot 
and cold days, network congestions and trip of generation units etc. 
load curves profile does not change so much. Statistical point of 
view, correlation constants for these five days are shown in Table 1: 

As can be seen in Table 1 that load curve correlations between 
these days are very high in positive direction and have average 0.96 
value. 

The same characteristic structure of load curve that is placed in 
Fig. 4 can be seen in hourly spot prices as well. The spot prices for 
electricity are generally low for the first 8h of the day. There are 
dramatical increase in prices after 7 a.m. while decrease in prices 
after 10 p.m. like the amount of electricity consumed by customers 
that has been previously shown in Fig. 4. 

In our case, because of above mentioned reasons, instead of 
choosing one spot price to represent risky asset in electricity mar- 
ket, it has been chosen 24h of day’s district spot prices to represent 
24 risky assets in electricity spot market. For analysis, hourly 
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Fig. 4. Daily load curve belongs to different weekdays between August of 2006 and 
December of 2009. 


system marginal prices which had been formed between August 
of 2006 and December of 2009 and hourly day-ahead spot prices 
which had been formed between December of 2009 and April of 
2011 are used as data sets. This time interval spreads 57 months 
and includes 1713 prices’ data for every hour of a day. That means 
41,112 prices’ data had been used in analysis. 


1,1 a2,1 24,1 
a4,2 a2,2 24,2 
f= : » h= ps a= : (8) 
41,1713 42,1713 24,1713 


Price vectors are formed as seen in Eq. (8). Where fn [n=1,2,..., 
24] represent price vector for each hour of day and an,m (n=1, 2, 
..., 24), (m=1, 2,..., 1713) represent hourly system marginal and 
day-ahead prices. 

In financial literature rate of return of an asset can be determined 
as the percentage change of investor’s wealth from the beginning 
to the end of a term. With electricity marketing point of view this 
approach is not true; we cannot economically store electricity or 
hold like other securities for a while to sell later and this product 
must be generated and consumed instantly. So it must be taken 
into consideration generation cost of electricity. So Min Liu and 
Felix Wu had determined rate of return for electricity spot market 
as; “rate of return=(spot prices — generation cost)/generation cost” 
[25]. Depends on your assumptions this generation cost param- 
eter can be taken as dynamic or passive. In dynamic approach 
CŒ = Sia (a + bp, + cpk?) tw ) power of generation unit P, and wt 
fuel prices are used as a parameter [24]. Besides C* = Sk + Dy py equa- 
tion can also be used with constant cost S;, variable costs D, and 
power of unit px to represent cost term of calculation. In passive 
approach like in this study it can be taken into consideration real 
average total electricity generation costs for a generation unit as a 
constant term. So that all price vectors in Eq. (8) are transferred to 
the rate of return vectors in Eq. (10) by applying the calculations in 


Eq. (9). 


eee (n=1,2,...,24), (m=1,2,...,1713) (9) 
11,1 12,1 T24,1 
11,2 12,2 T24,2 

n=| oi |, m= pee T=] | (10) 
11,1713 T2,1713 T24,1713 


where C’ represents marginal electricity generation costs and “*” 
index on it represents generation type (H: Hydro Power, T: Thermal 
Power with lignite and N: Natural Gas Combined Cycle Table 2). In 
this study EUAS’s real average electricity generation cost values has 
been used as a generation cost data set. rnm (n=1, 2,...,24),(m=1, 
2,..., 1713) represent rate of returns and rp is rate of return vector. 
After determination of these vectors, the average rate of return and 
variances of each set are easily calculated by using Eqs. (11) and 
(12): 


1713 


Š 1 
™ = 7993 Sora (11) 
m=1 
1 1713 
_ a _ #4 
On = | 7943 2 nm — Fn) (12) 
m=1 
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Table 1 

Correlation matrix for five independent days. 
Dates 21.09.06 20.03.07 08.06.07 07.01.08 18.11.09 
21.09.06 1 0.9514011 0.9703129 0.9706092 0.9684779 
20.03.07 1 0.9251289 0.9361949 0.9470435 
08.06.07 1 0.9605889 0.9486646 
07.01.08 1 0.9853352 
18.11.09 1 


4.2. Construction of three base models 


Before developing base form of Markowitz mean-variance opti- 
mization model, constitution of Variance-Covariance Matrix is 
necessary. By using Eq. (13), all elements of Variance—Covariance 
Matrix in Eq. (14) are determined as follows: 


1713 


1 - X 
Oxy = ig em — Fx)(Ty,m — Ty) (13) 
m=1 
o? 01,2 01,24 
o1 of m 
024,1 Oi 


These results will be used three main model approach. 


4.2.1. Standard mean-variance optimization model 

In this case electricity generators have only opportunity to 
sell their generations to spot electricity markets. That means it 
is constructed and used optimization model with 24 risky selling 
alternatives. It is called as standard mean-variance optimization 
approach and the mathematical programming formulations are 
formed as follows: 


24 24 
Min. (o2) = 5 XiX oij (15) 
i=1 j=1 
s.t. 
24 
tS SoXnin =Te (16) 
n=1 
24 
XXn 2 (17) 
n=1 
0O<Xı<1, [n=1,2,...,24] 


In above mentioned equation set electricity generators can sell, 
whatever amount of electricity they want to sell, in every hour of 
day without limitation. Firstly Eqs. (15)-(17) are used for formation 
of possible efficient portfolio alternatives and drawing of efficient 
frontier. Then, in second step is here finding the optimal portfolio 
with in efficient portfolios alternatives with respect to electricity 


Table 2 
EUAS average electricity generation costs according to generation type 2010 [26]. 


Indifferences Curve Efficient Frontier 
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Fig. 5. Portfolio optimization results for risky assets. 


generators’ risk aversion characteristic. To find the optimal portfo- 
lio, utility function (U) as mentioned before in Eq. (7) is used and 
optimization equation set is formed as follows: 


Maxy, U = E(rp) — Z402 (18) 


s.t. 
24 
SoXn = 1 O<Xn<1, YXne[n=1,2,...,24] 
n=1 


where 


E(tp) = X Xnřn (19) 


o2 = X XXa (20) 


To represent average risk-aversion level, “A” is taken as “3” 
[2,25]. Some researchers in financial literature has been tried to 
guess the value of “A” for average risk aversion. In spite of there is 
no strict value determined for “A” it is accepted has a value between 
2 and 4 for average risk aversion. After solution of above determined 
equations set it is reached optimal weight percentage of risky assets 
in the optimum portfolio for an electricity generator that is only 
taking into consideration of 24 risky alternatives in spot market to 
sell electricity. To reach the minimum global variance portfolio it 
is necessary the solution of Eq. (21). After solutions of all equations 
in this case it is reached result graphic as shown in Fig. 5: 


. 1, 2 
lim { Maxx, U = E(rp) - aR) (21) 


Type of power plant Generation (kWh) 


Industrial cost (cent/kWh) Commercial cost (cent/kWh) 


Natural gas combined cycle 
Lignite fired thermal power 
Hydraulic power plants 


20,107,683,691 
26,429,805,243 
39,914,183,522 


10.15 10.31 
7.61 8.50 
0.75 0.99 


Electricity Generation Co. Inc. (EUAS) 2010 real average electricity generation cost values. EUAS is a state owned company with the biggest market share (approximately % 


45 in generation and %50 in installed capacity) in Turkish electricity market. 


Commercial cost includes industrial cost and all other costs like administrative expenses, marketing, operating and maintenance costs. 
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4.2.2. Mean-variance optimization model with upper investment 
constraints 

Unlike the standard Markowitz mean-variance approximation, 
in this model, electricity generators cannot sell all energy in one of 
spot market hour they want. Instead of this, there are some con- 
straints for electricity generators for selling energy. This model 
basically is more realistic than standard model. In real electric- 
ity market environment generators (strategies can be change with 
respect to type of generators; hydraulic, coal, natural gas, nuclear, 
etc.) generally prefer less risky and stable conditions. Selling all 
of their energy to spot market without limitation is not a realistic 
case. So that it is constructed and used optimization model with 24 
risky selling alternatives with upper constraints for selling energy 
to each hours and equation sets are formed as follows: 


24 24 
Min. (o2) = DRG (22) 
i=1 j=l 
Si 
24 
Tp = XO Xnřn = Te (23) 
n=1 
24 
X X= (24) 
n=1 
Xn <On, [n=1,2,..., 24] (25) 
0<Xn, [n=1,2,...,24] 


In above equation sets it is used same standard equations with 
standard model (Eqs. (15)-(17)) we talked about in case a but there 
are limitations Xn < n for selling electricity in spot market hours in 
this model. These limitations can be gotten different values also but 
in our model all of them are taken same and equal to 0p. After the 
solution of Eqs. (22)-(25) efficient portfolios and efficient frontier 
are obtained like in previous case. Second step is here finding opti- 
mal portfolio with in efficient portfolios alternatives with respect 
to electricity generators’ risk aversion characteristic. To find the 
optimal portfolio, utility function as mentioned before in Eq. (7) is 
used and optimization equation set is given by as follows: 


1 


Max x,, U = E(Tp) = 3^0 (26) 
s.t. 
24 
Some 
n=1 
Xn <On, [n=1,2,..., 24] 


0<Xn, VXne[n=1,2,...,24] 


where 
24 
E(rp) = SOXnin (27) 
n=1 
24 24 
o2 = X XXa (28) 
i=1 j=1 


To represent average risk-aversion level “A” is taken as “3”. After 
solution of above determined equation sets it is reached optimal 
weight percentages of risky assets in the optimum portfolio for an 
electricity generator that is taking into consideration of 24 risky 


alternatives in spot market to sell electricity with 6, limitation 
for each hour. Reaching to minimum global variance portfolio it 
is also necessary solution of Eq. (21). After solutions of all equa- 
tions placed above it is reached results graphic as shown in Fig. 5 
same as previous case. 


4.2.3. Mean-variance optimization with one risk-free asset 

This case includes 24 risky assets and additional one risk-free 
asset for electricity generators. Bilateral contracts under guarantee 
of clearing authorities assumed as risk-free assets with their fixed 
prices for a fixed period of time. There are no limitations for elec- 
tricity generators to sell their electricity in risky spot market hours 
or via bilateral contracts. That means it is constructed and used 
optimization model with one risk-free and 24 risky selling alterna- 
tive. It is called as mean-variance optimization with one risk-free 
asset approach and related equation sets are formed as follows: 


24 24 
Min.(o2) = SOS CX Xjo% (29) 
i=1 j=1 
s.t. 
24 
ee SoXain +Xytp = Te (30) 
n=1 
24 
XXn +X =1 (31) 
n=1 
0< Xr, <1 


0<Xn<1, VXne[n=1,2,...,24] 


By solution of these equation sets for every expected rate of 
return it is obtained efficient portfolios and efficient frontier. Eq. 
(31) shows total weight percentage of portfolio assets including 
risk-free asset and this value cannot exceed 1. There are no limita- 
tions for electricity generators like in case a but in this case extra 
one risk free asset (bilateral contract option) for electricity gener- 
ators. To prevent from risk, more risk averse generators can prefer 
selling their generations into bilateral market instead of spot mar- 
ket. To find the optimal selling portfolio for electricity generators 
with average risk-aversion, value of “A” is taken 3 again and Eq. (7) 
is used for utility function (U) and the optimization constraints as 
follows: 


Maxx, xy U = E(re) — sho? (32) 


s.t. 
24 
XXn +X, =1 
n=1 


0<Xn<1, VXne[n=1,2,...,24] 


0< Xr, <1 

where 

E(rc) = Xrerf + (1 — Xr )E(rp) (33) 

o2 = (1 — Xy 03 (34) 
24 

Erp) = X Xnřn (35) 
n=1 
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Table 3 
Risky assets for hydraulic power plants. 


H Expected rate of return fn (%) Standard deviation on (%) H Expected rate of return fn (%) Standard deviation on (%) 
1 891.68 239.52 13 982.25 275.93 
2 768.03 300.31 14 997.67 277.90 
3 624.91 328.41 15 1019.46 309.23 
4 518.24 324.42 16 995.17 302.51 
5 472.82 308.34 17 987.67 288.57 
6 451.67 291.21 18 971.91 357.12 
7 412.64 302.78 19 932.62 315.31 
8 598.28 366.44 20 917.25 319.01 
9 865.04 342.74 21 929.85 310.05 
10 962.99 291.16 22 912.71 322.65 
11 1011.41 262.98 23 970.53 195.29 
12 1044.88 275.02 24 930.89 220. 84 
Table 4 
Risky assets for lignite fired thermal power plants. 
H Expected rate of return rp (%) Standard deviation on (%) H Expected rate of return rp (%) Standard deviation on (%) 
1 15.50 27.90 13 26.05 32.14 
2 1.10 34.98 14 27.85 32.37 
3 —15.57 38.25 15 30.38 36.02 
4 —27.99 37.79 16 27.55 35.23 
5 —33.28 35.91 17 26.68 33.61 
6 —35.75 33.92 18 24.85 41.59 
7 —40.29 35.26 19 20.27 36.72 
8 —18.67 42.68 20 18.48 37.16 
9 12.40 39.92 21 19.95 36.11 
10 23.81 33.91 22 17.95 37.58 
11 29.45 30.63 23 24.69 22.75 
12 33.34 32.03 24 20.07 25.72 
Table 5 
Risky assets for natural gas combined power plants. 
H Expected rate of return rp (%) Standard deviation on (%) H Expected rate of return rp (%) Standard deviation on (%) 
1 —4.48 23.00 13 3.92 26.50 
2 —16.65 28.84 14 5.40 26.69 
3 —30.39 31.54 15 7.49 29.69 
4 —40.63 31.15 16 5.16 29.05 
5 —45.00 29.61 17 4.44 27.71 
6 —47.03 27.96 18 2.93 34.29 
7 —50.77 29.07 19 —0.84 30.28 
8 —32.95 35.19 20 —2.32 30.63 
9 —7.33 32.91 21 —1.11 29.77 
10 2.07 27.96 22 —2.76 30.98 
11 6.72 25.25 23 2.80 18.75 
12 9.93 26.41 24 —1.01 21.21 
24 24 with CY = 10.31 cent/kWh. They have been calculated means and 
op = 5 XiXjoij (36) variances of each hour for three special cases and shown in the 
i=1 j=l following tables (Tables 3-5). 


where Xp is the weight percentage for risk-free asset (bilateral con- 
tract), E(rc) is the expected rate of return for mixed portfolio which 
constitutes of risky and risk-free assets, and o2 is the variance of 
mixed portfolio. 

By solution of above equation sets it can be obtained weight 
percentages of optimal portfolio in a risky and risk-free market 
environment. To represent average risk-aversion level, “A” is taken 
as 3 again. Reaching to minimum global variance portfolio it is also 
necessary solution of Eq. (21). After solutions of all Eqs. (29)-(36) it 
is reached results graphic as shown in Fig. 5 same as previous cases. 


5. Results 


All cases which have been determined in previous sections 
are applied to three specific examples; Hydraulic Power Plants 
with C" =0.99 cent/kWh, Lignite Fired Thermal Power Plants with 
C7 =8.50 cent/kWh, and Natural Gas Combined Cycle Power Plants 


In all cases while finding optimum portfolios, “A” is taken 3 for 
average risk aversion and 1000 for global minimum variance port- 
folios. Actually with respect to Eq. (16) “A” must be taken infinity 
to find global variance portfolios but in practice this is not appli- 
cable situation and it has been seen some error propagations had 
occurred when taken “A” with high values. So “A” has been taken 
1000 to find global variance portfolios. 


5.1. Efficient frontier and optimum portfolio for hydraulic power 
plants 


Case a. Standard mean-variance optimization: there are 24 risky 
assets for electricity generators and no limitation for selling of elec- 
tricity in spot market’s hours. “OP” represents optimum portfolio 
and “GMVP’” represents global minimum variance portfolio in Fig. 6 
and next graphics. 

Case b. Mean-variance optimization model with upper invest- 
ment constraints: there are 24 risky assets for electricity generators 
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Fig. 6. Portfolio optimization for hydraulic power plants with 24 risky assets. 
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Fig. 7. Portfolio optimization for hydraulic power plants with 24 risky assets and 
upper investment constraints. 


but there are also upper limit determined for electricity trading for 
each hour. When this limit is taken as 12.63% for each hour it has 
been reached results in Fig. 7. 

Case c. Mean-variance optimization model with one risk-free 
asset: in this case beside the 24 risky assets there is also one risk- 
free asset (bilateral contract). Rate of return for bilateral contract is 
assumed as 750% (8.415 cent/kWh, this is because of low generation 
costs of hydraulic power plants) (Fig. 8). 


5.2. Efficient frontier and optimum portfolio for lignite fired 
thermal power plants 


Case a. Standard mean-variance optimization: there are 24 risky 
assets for electricity generators and no limitation for selling of elec- 
tricity in spot market’s hours (Fig. 9). 
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Fig. 8. Portfolio optimization for hydraulic power plants with 24 risky and one risk- 


free asset. 
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Fig. 9. Portfolio optimization for lignite fired thermal power plants with 24 risky 
assets. 
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Fig. 10. Portfolio optimization for lignite fired thermal power plants with 24 risky 
assets and upper investment constraints. 


Case b. Mean-variance optimization model with upper invest- 
ment constraints: There are 24 risky assets for electricity generators 
but there is also upper limit (12.63%) determined for electricity 
trading for each hour (Fig. 10). 

Case c. Mean-variance optimization model with one risk-free 
asset: In this case beside the 24 risky assets there is also one risk- 
free asset (bilateral contract). Rate of return for bilateral contract 
is assumed as 8% (9.18 cent/kWh, because of high generation costs) 
(Fig. 11). 
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Fig. 11. Portfolio optimization for lignite fired thermal power plants with 24 risky 
and one risk-free asset. 
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Fig. 12. Portfolio optimization for natural gas combined cycle power plants with 24 
risky assets. 
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Fig. 13. Portfolio optimization for natural gas combined cycle power plants with 24 
risky assets and upper investment constraints. 


5.3. Efficient frontier and optimum portfolio for natural gas 
combined cycle power plants 


Case a. Standard mean-variance optimization: there are 24 risky 
assets for electricity generators and no limitation for selling of elec- 
tricity in spot market’s hours (Fig. 12). 

Case b. Mean-variance optimization model with upper invest- 
ment constraints: there are 24 risky assets for electricity generators 
but there are also upper limits (12.63%) determined for electricity 
trading for each hour (Fig. 13). 

Case c. Mean-variance optimization model with one risk-free 
asset: In this case beside the 24 risky assets there is also one risk- 
free asset (bilateral contract). Rate of return for bilateral contract is 
determined as 4% (10.7224 cent/kWh, because of very high gener- 
ation costs) (Fig. 14). 


5.4. Analysis of results 


In case a, similar results have been obtained for all electricity 
generators with different weight percentages for same hours. 


Table 6 


10.00% 


0.00% 
% 
ie 0.3 9035 04 
= -10.00% 
E 
s 
2 
© -20.00% 
E 
[e] 
2 9, 
-30.00% è 
a + 
-40.00% é 


-50.00% 


Risk (0) 


Fig. 14. Portfolio optimization for natural gas combined cycle power plants with 24 
risky and one risk-free asset. 


As can be seen in Table 6, weight percentages of hour-12 in opti- 
mum portfolios for electricity generators increase with increasing 
electricity generation costs. Optimum portfolios for lignite and nat- 
ural gas fired electricity generators are same but have different 
weight percentages for hour-12 and hour-23. In hydraulic there 
is only one option with 100% of hour-23. Under normal conditions 
for all three types of electricity generators, most profitable spot 
market hour (highest average spot prices) is hour-12 but weight 
percentage of this hour for electricity generator that has hydraulic 
power plant is only 0%. The reasons of results must be evaluated 
according to risk and return concept. In our application standard 
Markowitz mean-variance optimization model in case a has been 
successfully applied to all electricity generators and has generated 
meaningful results. Assuming that expected rate of returns have 
normal distribution, z tests have been applied to results that are 
placed in Table 6. Possibilities of getting positive rate of return for 
optimum portfolio preferences of each electricity generators have 
been obtained. These values are approximately 100% for hydraulic 
power plants, 87.9% for lignite fired thermal power plants and 63.8% 
for natural gas combined cycle power plants. The advantage of low 
generation cost has been easily seen from the results. 

In case b, different results have been obtained for all electricity 
generators (Table 7). 

In this model called as case b, upper investment constraints have 
been applied (selling limit for spot market hours) to all electricity 
generators. Results that have been obtained for trading hours 11, 
12, 14, 23, 24 are common with 12.63% weight percentages. This 
case is more realistic than previous case a. Because in real market 
environments it cannot be sold all available energy to one market 
hour, there must be some constraints for selling electricity in spot 
market to real illustration of electricity market environment. It is 
also necessary to be aware of characteristics of different types of 
power plants. For example, a hydraulic power plant has a chance to 
wait for a while before selling electricity because it has a reservoir 
holding water in it. Same case is not valid for thermal power plants. 
Start-up and shut-down processes are very expensive and it needs 
time for smooth cooling and warming to prevent thermal power 
plant from thermal shocks. An average start-up time required for 
a lignite power plant is between 8 and 12h for hot start-up and 


Optimum portfolio results of standard mean-variance optimization model according to types of electricity generators. 


Optimum portfolio Hydraulic Lignite fired Natural gas 
Electricity generation cost 0.99 cent/kWh 8.50 cent/kWh 10.31 cent/kWh 
Hour 12 - 58.08% 72.49% 

Hour 23 100% 41.92% 27.51% 
Expected rate of return 970.53% 29.71% 7.97% 

Standard deviation 1.952 0.266 0.233 
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Table 7 
Optimum portfolio results of mean-variance optimization model with upper investment constraints according to types of electricity generators. 
Optimum Portfolio Hydraulic Lignite Fired Natural Gas 
Electricity generation cost 0.99 cent/kWh 8.50 cent/kWh 10.31 cent/kWh 
Hour 1 12.63% - - 
Hour 2 0.16% - - 
Hour 6 3.99% - - 
Hour 7 10.14% - - 
Hour 10 1.97% - - 
Hour 11 12.63% 12.63% 12.63% 
Hour 12 12.63% 12.63% 12.63% 
Hour 13 6.94% 11.59% 11.59% 
Hour 14 12.63% 12.63% 12.63% 
Hour 15 - 12.63% 12.63% 
Hour 17 1.02% 12.63% 12.63% 
Hour 23 12.63% 12.63% 12.63% 
Hour 24 12.63% 12.63% 12.63% 
Expected rate of return 896.78% 27.33% 4.97% 
Standard deviation 2.176 0.287 0.236 
Table 8 
Optimum portfolio results of mean-variance optimization model with one risk-free asset according to types of electricity generators. 
Optimum portfolio Hydraulic Lignite fired Natural gas 
Electricity generation cost 0.99 cent/kWh 8.50 cent/kWh 10.31 cent/kWh 
Bilateral contract price 8.42 cent/kWh 9.18 cent/kWh 10.72 cent/kWh 
Hour 12 6.08% 58.08% 28.38% 
Hour 23 12.73% 41.92% - 
Bilateral contract 81.20% - 71.62% 
Expected rate of return 795.99% 29.71% 5.68% 
Standard deviation 0.392 0.266 0.075 


maybe larger than 1 day for cold start-up. Natural Gas Combined 
Cycle Power Plants need less time with respect to lignite fired ther- 
mal power plants but still it is expensive process for natural gas 
combined power plant, too. 

In this application mean-variance optimization model with 
upper investment constraints (case b) has been successfully applied 
to all electricity generators and has been generated meaningful 
results. Assuming that expected rate of returns have normal distri- 
bution, ztests have been applied to results that are placed in Table 7. 
Possibilities of getting positive rate of return for optimum portfo- 
lio preferences of each electricity generators have been obtained. 
These values are approximately 100% for hydraulic power plants, 
82.89% for lignite fired thermal power plants and 58.32% for natural 
gas combined cycle power plants. 

In case c; in the final case to reach more realistic results, one 
risk-free asset has been taken into consideration besides 24 risky 
assets. 

Bilateral contract prices have been generated according to 
reasonable profit margin. These prices have been assumed as 
8.42 cent/kWh for hydraulic power plants, 9.18 cent/kWh for lig- 
nite fired thermal power plants and 10.72 cent/kWh for natural 
gas combined power plants. Case c has also been applied to all 
three electricity generators. Results placed above in Table 8 have 
been obtained. Results have similarities with case a except bilateral 
contract weight in optimum portfolio but for lignite fired ther- 
mal power plants same results have been obtained in case a and 
case c. Determination of profit margin is very important part of 
this application case. Because profit margin directly affects the 
weight percentages result of optimum portfolios. High electricity 
generation costs and low profit margins cause to reduction of spot 
market share in optimum portfolio to reach acceptable level of rate 
of return and risk. Hydraulic power plants with their low genera- 
tion cost and generation program which strongly depends on water 
regime do not prefer bilateral contract so much. Average capac- 
ity factor for hydraulic power plants are assumed as 40% and can 
change from country to country because of geographical reasons. 
But the situation is different for lignite power plants and natural gas 


combined cycle power plants. To reduce non-generation risk they 
need bilateral contracts in their selling strategies and in our analy- 
sis also shows that to reach acceptable rate of return and risk levels, 
bilateral contract’s share must be increased in optimum portfolios. 


6. Conclusions 


This study provides an investigation on the significance of the 
mean-variance optimization technique for electricity markets in 
Turkey. It formulates the general portfolio optimization problem 
for Turkish day-ahead electricity market. Three cases were estab- 
lished with different approaches like in financial literature. It was 
shown that the solutions of portfolio optimization problems for “n 
risky assets”, “n risky assets with upper investments constraints” 
and “n risky assets with one risk-free asset” for hydraulic, lignite 
fired thermal, natural gas combined cycle power plants of EUAS 
which is the biggest electricity generation company of Turkey. 

As a result, mean-variance optimization methodology with 
three different cases can be applied to energy allocation between 
spot market and bilateral contracts in a market where there is no 
transmission congestion in the system. In Turkey, there are no local, 
zonal or nodal pricing systems so that spot market price is used as 
a one price everywhere in market. Therefore, it has not been taken 
into consideration risks because of congestion charge. Turkish bal- 
anced market and day-ahead planning market data are used to 
illustrate the day-ahead Turkish electricity market to apply mean- 
variance optimization. Price occurred every hour is thought as one 
risky asset so that it can be determined 24 risky selling alterna- 
tives for electricity generators and this is a new assumption such 
a kind of study. Electricity generators depends on the constraints 
(production program, maintenance program, fuel prices, network 
congestions, etc.) of course and by taking into consideration risk 
and rate of return preferences can prefer to sell their electrical 
energies in different hours or via bilateral contracts. 

The methodology showed here includes mean-variance opti- 
mization with three cases for three different electricity generators. 
This methodology can also be applied to more sophisticated and 
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detailed problems. Fuel prices, transmission congestions, main- 
tenance program, trips, transaction costs and financial markets 
(option, swap, forward and futures) are be used as a factor for 
improvement of this approach. Other risk management techniques 
hedging, risk measurement and asset valuation can also be applied 
to Turkish electricity market. These are determined as prospec- 
tive studies. Consequently, it has been proved the applicability of 
mean-variance optimization, which is famous in finance literature, 
to Turkish electricity market. 
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